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I. Estabrook-Wahlquist Prolongations and Zero-Curvature Requirements 

Motivated by a desire to find new solutions of physically-interesting partial differ- 
ential equations, we think of a k-th order pde as a variety, Y, of a finite jet bundle, 
(^jVf, N), with M the independent- and N the dependent- variables for the pde. From 
•/^ ■ this geometric approach, we can look for point symmetries, contact symmetries, gener- 
Q\ • alized symmetries, or even non-local symmetries, where the system is prolonged further, 
to a fiber space over J°°, with fibers W, where vertical fiows map solution spaces of one 
Qj [ pde into another, satisfied by the additional dependent variables, w"^, that coordinatize 
Q I the fibers. The compatibility conditions for such fiows to exist are "zero- curvature con- 
^ ' ditions. " Solutions of these conditions may be found using the tangent structure or the 

co-tangent structure, over J°° x W. We describe both, but follow the approach via dif- 
! ferential forms, following Cartan,^ , Estabrook and Wahlquist,^ and Pirani,"^ believing 
that it provides better guides for the intuition, for complicated (systems of) pde's. 

For a vector-field presentation, we choose a commuting^ basis, {ca}, for tangent 
vectors over M, and lift them to the total derivative operators. Da, over J°°. Provided 
the system of pde's is involutive, they will still commute when restricted to the variety 
described by the pde's, this restriction being denoted by Da- The further prolon- 
gation into the fibers W requires the addition of some vector fields vertical with respect 
O ■ to the fibers, i.e., = X^id/dw"^), with the functions of both the jet variables 
and the {w^}. Requiring that the Da + X„ commute, when restricted to Y^^^ x W, for 
^ ; some £, ensures that the can act as pseudopotentials for that pde:^'^ 
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[Da + X^,D, + X,]^ ={Da{X^)-D,{X^)}-^ + [X^,X,] . (1.1 



The general solution for the describes all possible fiber spaces, or coverings,^ for this 
pde. As the construction gives the Xa the "form" of a connection, it is reasonable to refer 
to these equations as "zero-curvature" requirements; it is, however, a generalization of 
the usual approach, ^'^ since the Xa's are still only elements of an abstract Lie algebra 
of vector fields, with neither coordinates, nor even their number yet determined. 

As an identity in the jet coordinates, Eqs. (1.1) determine several independent 
equations. Their solution describes the X^ as linear combinations of vector fields W^, 
with coefficients depending on coordinates for Y^^\ with the w'^-dependence encoded 
within a set of commutation relations among the {Wq,}, as vector fields within the 
algebra of vector fields over W. The smallest subalgebra that faithfully reproduces the 



linear independence, and the values, of those commutators is the general solution to the 
covering problem; we believe it is a universal object for the given pde and others related 
to it, so that it may be used to characterize related classes of pde's.^'^'^° 

The isolation and identification of such algebras is an important part of the pro- 
cess of determining and understanding all the solutions of nonlinear pde's. Vector- 
field realizations will generate Backlund transformations, inverse scattering problems, 
etc.-'^^'-'^-'^ Faithful realizations will usually involve infinitely many pseudo-potentials, 
making their identification somewhat difficult, and the first researchers did not con- 
sider the entire infinite-dimensional algebras. However, beginning with the work by 
van Eck,-*^^ and Estabrook,-*^^ on identification of the universal algebra for the KdV 
equation, the search for the infinite-dimensional algebras involved has been extended 
considerably by Hoenselaers and co-workers^^'^'^'-'^^ , by Omote^^ , and by the group 
at Twente, who seem to have made this a studied art-form. The dual approach, via 
differential forms, created by Estabrook and Wahlquist and built on the ideas of Cartan, 
begins with the 'contact module,' f]'^(M, N) C [j'^(M, N)]^, generated by the following 
set of 1-forms: 

Q'^{M,N):\ ... \={e^\\a\ = 0,...,k-l} 

\ "aia2---ak-i "'^ai02...afc_i ^aia^-.-ak-iak""^ J 

(1.2) 

where the summation convention is being used, and a choice for a local coordinate chart 
is {x", 2^, • • • 5 -^ai...afe with standing for any of these (jet) coordinates 

except the independent variables themselves, x"". The contact module 'remembers' the 
relation the coordinates of the jet bundle would have when they are pulled back by 
the lift of a function over M: u : U C M ^ N ^ (/w)*(0'^) = 0. The ideal, X, 
is the differential closure of the puUback of the contact module to Y, and constitutes 
the Cartan description of the original pde. For 2 independent variables, {x,y}, the EW 
procedure first chooses a proper, closed subideal, /C C X, generated by a set of 2-forms, 
{a^}, that still is effective at describing the given pde.^®'^^ The new variables {tu^} are 
adjoined by appending to /C contact forms, {<^^}, for each of these pseudopotentials, 
and maintaining the ideal closed: 



u^ = - dw^ + F^dx + G^dy, 
dF^ A dx+dG^ Ady = fi + r/^^ A 



N . (1.3) 



To show "equivalence" with the zero- curvature equations, we first consider all of 
J^'^\ i.e., without a pde, and then "restrict" to Y. For simplicity considering quasi-linear 
pde's, we first select £, generated by the wedge-products of all 1-forms in with the 
dx"". For < < k — 1, C, contains exactly one copy of each of dz^. Labelling its 
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coefficient, in Eqs. (1.1), by {f^)%, we have 

F% = {fK, , Gl^ = (/<^ , (1-4) 

for each jet coordinate z^, with no repetitions. Writing {r]^B)z^ ^r the components of 
the 1-forms ij^b, Eqs. (1.4) also gives us 

= (?7^b)u,c dw'^ A dw^ , = [r]^ b) z'^dz'^ A dw^ , 
F^sdw^ ^dx= -{r]'^B)xdx ^dw^ , C^^sdw^ A dy = -{r]'^ B)ydy A dw^ , (1.5) 
=^ ri^B^F^sdx+G%sdy . 

The only remaining part of Eqs. (1.4) are the coefficients of the basis 2-form dx A dy: 

-F^y + Gi = -F^ + Ft. + J2 {-''oyF% + z^^,G% } . (1.6) 

kl=o 

Introducing vertical vector fields, F = {F'^)d / dw"^ and G = {G^)d/dw^, so that the 
first two terms on the right hand side are the components of a commutator, this becomes 

k-i 

[D, + F,Dy + G]^ [d, + F,dy + G]+Y, l-^^yF,.- + ^^.G,,.} = . (1.6') 

kl=o 

This has the same form as Eqs. (1.1), except that we must still effect the restriction 
to some variety Y. The resulting EW ideal, /C, will be defined over 
j(fe-i)^ so that it reduces the problem in an important way; F and G will depend on 
several fewer variables — only the coordinates for Y^'^~^\ which we select by choosing 
"co-coordinates" for Y, i.e., those jet coordinates the pde's will be used to eliminate, 
as a method of (locally) defining Y C j(^).20 if the restriction of £ to y removes all 
the k-th level coordinates, then it can be taken as /C. Otherwise, the remaining highest 
derivatives must still be removed from the system, which is always possible, although the 
methods depend on the particular pde. For a simple, quasi-linear evolution equation, 
we choose Zy — H Z(^k) + K as our co-coordinate, where H and K are functions over 
Y^'^~^\ Restriction to Y then causes Z(^k) to appear twice: in (dz — Zy dy) A dx and also 
((i^(fe_i) — dx) A dy. The following replacement process, followed by dropping the 
second 2-form above, reduces £. to {yC^"^)}*, as desired: 

(dz — Zy dy) Adx ^ {^dz — {H Z(k) -\- K) dy^ A dx 
=dz Adx + H dz(^k_i) Ady — K dy Adx mod {dz(j^_i-) — 2;^^) dx) A dy . 

For an evolution equation, the above process is unique, and the EW process gives 
exactly the same results as that using vector fields; however, in general the situation 
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is different. As a second example, we consider a pde that defines Y via the equation 
Zyy = H Z(^i^^ + J Z(^k-i),y + K, where the integers in parentheses indicate the number 
of a;-derivatives, and H, J, and K are defined over Y^^~^\ The earher replacement 
process now has more than one path allowed: 

{dzy — Zyy c/j/) A dx — > dzy A dx — K di/ A dx + H dz(^k-i) A dy 

dz^k-2),y ^ dy ' Option 1 

^ I — <^^(7,_i) A dx , Option 2 (equiv. to 1) i^-^) 

\^'^{k -1) ^dy — a dx) , Option 3, symmetric, 

but inequivalent. 

The two inequivalent paths each generate acceptable EW ideals, but lead to distinct 
prolongation structures. A third example pde with inequivalent prolongation structures 
is the sine-Gordon equation, as described below. 

We describe the generators in /C as a set of contact 2-forms for only those coordi- 
nates needed for y('^~'^\ and a (set of) "dynamical 2-forms" for each pde in the system. 
The ideal /C provides us a geometrically-motivated structure for knowing on which jet- 
variables we need no dependence; as well, the Lagrange multipliers, f^rj expressed in 
terms of derivatives of the F"^ and G^, as in Eqs. (1.4), tell us on which of the coor- 
dinates they must depend. (The curvature should vanish only when it is restricted to 
y, so that the are truly pseudopotentials for the pde; within the functional form 
of the curvature, the f^r multiply the 2-form expression of the pde. The remaining 
information is then the commutator equation, Eq. (1.6'). 

II. Uxy = f(x, y;u) : The sine-Gordon and Robinson- Trautman equations 

Using the generalized form of the equation, with any choice of /, we begin with the 
ideal, £, as just described. Using the pde to replace u^y within £ leaves us with the 
largest possible EW ideal, with 4 generators. This ideal generates exactly the vector-field 
commutator equations,^ and is too large to allow us to solve the resulting equations. 
However, Pirani,^ has shown the existence of two distinct, useful sub-ideals: 

{( {du — p dx) Ady ^ 
(du — pdx) A dy \ 
L IC^: { {du-qdy)Adx ) . (2.1) 

{dp - fdx) Adx \ 

\^dp Adx - dqAdy + 2fdxAdy ) 

For the case f = f{u), only, such as the sine-Gordon equation, we have 

[F , G] = G„ + g F„ + f{u) (F^ - G,) , 
with F = F(w,p;m;^) , G = G(w, zw^) , (2.2) 
and /Ci : F„ = = Gg , IC2 : F^ + Gg = . 
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These equations are of quite a different character than those for evolution equations. For 
an evolution equation, the constraints resolve the dependence over Y by the solution 
of algebraic equations; however, for higher order equations, vector-field-valued pde's 
must be solved, substantially increasing the difficulty of the problem.^-*^ Nonetheless, 
with relatively minor assumptions, the reductions generated by these smaller algebras 
allow us to resolve these equations as the "flow" of one vector field along the direction 
described by another, and to express the solutions explicitly in terms of the adjoint op- 
eration of one field upon another. We describe the general solutions of each, labelling 
the resultant algebras by 2li and 2t2, and the prolongations {F, G} within them as 
{^i, Qi \ i = 2}. We also show that 2I2 is gauge equivalent to a subalgebra of and 
that a subalgebra of that is homomorphic to A\^\ Important earlier work on infinite 
versions of these algebras was done by Hoenselaers,^*^'^^ by Leznov and Saveliev,^ by 
Dodd and Gibbon^"^ for Ai, and by Shadwick,^^ for A2- 

With details in Ref. 22, the solution for 2li may be found by first defining Z = 
Qi + {Gi)uu, which requires [JF, Z] = — pZ„. Expanding JF about the origin and writing 
F„ as the coefficient of p'^/nl, the fiow equations tell us that 

2 = e-«(adF,)Q^ , [F, , (adFirQo] = 0, , n ^ 1 , m = 0, 1, 2, . . . , 

(2.3) 

where Qo G is a "constant" of the integration. Continuing in this mode, the general 
solution is 



^i-Fo+pFi= / d.(p-.)e^(^d«^)Ko = J]/— -(adGifKo, 

PU AT? °° 

^i-GocosM-Gisinw= / dwsm{u - w)e-'" ^^^^^^Qo = V (-cosw)^-"-2^(adFi)''Qo 

•^0 n=0 

and requirements on (only) some of the commutators: 

Sq] = ' tfo, gj = -fi : [fo: 3n] = , [fi, gj = g^ , [fi, gj = -g^ + q^ - ko , 

K„ = (adGi)"Ko , Qm = (-l)'"(adFi)"^Qo. 

Alternatively, the equations that define 2I2 first tell us that 

rj^2 = |pR + B , ^?2 = -|9R+C, R„ = , 

2t2 : <^ (2.4) 
[[iR,B]=B, , [|R,C]=-C, , [B,C]=R/(w). 

Integration of the differential equations gives two new, vertical vector fields such that 

B = e+^"«^^^))E, c = e-^'^«^^^»J . (2.5) 
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Defining iterated commutators, = (+adR)™'E, = (— adR)"J, and setting c„ 
as the coefficients multiplying tt^/n! in the power series for /(«), the content of the 
commutator equation becomes simply 

[Efc-m, Jm] = CfeR , ym30<m<k . (2.6) 



Even when divided by the countably infinite set of relations in Eqs. (2.6), the free 
algebra generated by {J, R, E}, is so far unidentified as an already-studied algebra. In 
fact, those relations do not appear to be consistent with the usual sorts of gradings, 
so that some distinct approach to infinite-dimensional algebras may be required. We 
mention two quite different avenues that can be followed at this point. If one no longer 
requires all the to be linearly independent, then the homomorphic mapping to A^^^ 
can be demonstrated. Alternatively, maintenance of linear independence seems to lead 
in the direction of Toeplitz algebras of operators over Banach spaces. 

Viewing F and G as the components of a Lie-algebra- valued connection 1-form, F, 
over the covering spaces, a gauge transformation generated by a vertical vector field, 
say S, would transform F by = e**^*^^*^ ^^^F — d{tS). Choosing S = — ^R, transforms 
2I2 into that part of 2li generated by setting Kq = 0. 

{R — ^ Fi , En — > Fn , Jn — > Gn , 
T r , o T 

dl —>■ — <J1 , . . . , Jj — > — tij-2 ■ 

Hoenselaers' homomorphism-'^^ of part of 2li into A^^^ can now be extended to 2I2: 

m-2 (1) 



E^=(adR)"'E ^ (adFi)"'Fo 




for m even, 
for m odd, 

m + n even. 



[E„.E„] ^[(adF.rF„,(adF.)"F„lv ((-D-^jf , ™ + „odd. ' 

Triple commutators of the £^,'3 among themselves will generate elements of A\^^ at the 
third level, etc. While interesting, this homomorphism loses much information contained 
within the larger algebra; for example, \l/'(E^+2) = ~*'(Em) ? '^'i^k+2) — — ^'(Jfc)' 
it eliminates any information carried by the generators and K^. 

A quite distinct approach would maintain linearly independent at least those of the 
generators that appear in J^2 and ^2, which leads us to consider that subalgebra spanned, 
as a vector space, on the countable list of generators {R, E^,J„ | m, n = 0,1,...}, 
therefore requiring the double, triple, etc. commutators to be linear combinations of 
these, such as [Em,,En] = {A^m)^ Ji- This is a quite distinct approach from that which 
led io A^^\ where only (R, Eq, E^, Jq, Ji}, of this original set, were linearly independent. 
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and it was their commutators that generated the higher levels of the Kac-Moody algebra. 
Solving the constraints on the coefficients determining the linear combinations leads to 
the creation^^ of a Banach algebra of Toeplitz operators made from countable sums 
of those coefficients that define the commutators. This mode of thinking causes us to 
re-describe the countable set of repeated commutators, E^, in terms of a function E(t), 
defined within a Banach spaces of functions on the circle, along with a set of integral 
equations, reminiscent of Weiner-Hopf equations, for functions defined on x S^. 
The resulting prolongation forms, J^2 and Q^, would then be expressed as integrals over 
that circle: 

We may now return to the case where / does depend on the independent variables; 
an important example for this case is the Robinson-Trautman equation for diverging, 
non-twisting, Petrov type III solutions of the vacuum Einstein field equations: 

RT equation of type III: Ua:y = l{x + y)e-2" . (2.9) 

Using the symmetric subideal, /C2, in Eqs. (2.1), the RT version of Eqs. (2.2) is 

F = pZ + B , G = -qZ + C , Z„ = , 
[Z,C] = -Cu + Zy , [Z,B] = +B^ + Z^ , (2.10) 
[B,C] = Bj,-C, + (a; + y)e-2'^Z . 

Since the sine-Gordon equation is well-studied, we were interested in "all" of its 
subtleties. On the other hand, no interesting solutions of the RT equation are known; 
therefore we will look, first, for as simple a solution as possible. Obviously the pro- 
longation structure must depend on {x,y}. We have considered two complementary 
cases in Ref. 26, integrating there these three-term pde's that generalize our earlier 
flow equations. We also show there that (at least) subalgebras for each case are gauge 
equivalent. Here, we only follow the simplest case, which assumes that = = Ga;, 
thereby requiring that Z be independent of both x and y, and reducing our pde's to 
simple flows, which integrate as before. The additional assumption that all commuta- 
tors of Z with their (respective) constants of integration are parallel gives us a system 
quite similar to the generators for sl(2), but with {x, y}-dependence: 

B(x,«)=e+"(^dZ)R(a;) , C(y, zx) = e'-^^^ Z) s(^) ^ 
[Z, S] = S , [Z, R] = -R , =^ [R{x),3{y)] = {x + y)Z . 
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As the independent variables occur linearly, the simplest, non-trivial solution is for 
R and S to be linear polynomials, which causes the entire EW prolongation algebra to 
have a contragredient form: 

R(a;) = -fi - xf2 , S{y) = +62 + yei , 

[Z,ei]=ei , [Z,fi] = -fi ,i = l,2 , (2-12) 
[e2,fi] =0= [ei,f2] , [ei,fi] =Z = [e2,f2] • 

The three Lagrange multipliers are now proportional to {R(a;), S(y), Z}. We must 
therefore determine a realization of the algebra defined by these 5 generators, which 
maintains these three linearly independent. At some level, this is a simple task, since 
this is in fact just the simplest contragredient algebra of infinite growth, usually referred 
to as -^^2,^^ when one identifies our Z with its generator h. The generic contragredient 
algebra^^ has the (standard) form 

[ej, fj-] = 5yhj , [hi,hj] = , [hi.Bj] ^ AijBj , [h^, f,-] = -A^f,- , (2.13) 

where the Aij are elements of a (generalized Cartan) matrix A. Our algebra K2 is just 
the algebra described by the matrix ^ = ( ]^ |) , divided by its center, h\ — h2. However, 
since no realizations of K2 have yet been discovered, this is not the end of the task. 

The infinite growth of K2 is of course the difficulty To see how it affects the 
problem directly, we note that the Jacobi identity requires [Z, {(adei)"^e2}] — {n -\- 
l){(adei)"e2}, and also [fa, {(adeO^+^ea}] = -|m(m + l){(adei)'"e2}. The first 
equality shows us the growth of the dimension of the space of z-th level commutators, 
unless the objects {(adei)'^e2} were to vanish from some value of n onward. The 
second equality tells us that this would cause a downward cascade, leaving us with zero 
values for our Lagrange multipliers. To give precise definitions, we define an (integer)- 
graded Lie algebra as one that can be presented as a direct sum of subspaces, which 
the Lie bracket operation "preserves"; i.e., for Q = ®Gi, we have [Qi , Qj] C QiJ^j. If 
is the dimension of Qj, then the (Gel'fand-Kirillov) growth,^^ r, of Q is defined 

as lim„_,.oo {log (i^/ log n}. For K2 one finds that dn grows like 2"', so that the resulting 
growth is infinite. More recently, Kirillov^^ has introduced the notion of algebras of 
intermediate growth, where \ogdn grows like n^, for some < 5 < 1, and has shown 
that Vect{R^) is an algebra of intermediate growth, with 5 = m/{m + 1). We conclude 
from this that K2, which has 5 — 1, will not have a realization within Vect{K^) for any 
finite value of m. 

Nonetheless, the next step in the process of finding new solutions is to determine 
explicit realizations of this algebra, use the variables in the carrier space as pseudopoten- 
tials, pick out a Backlund transformation, take the one existing solution, and begin to 
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generate new ones, as has been done many times before with many other pde's. Since 
this is indeed the minimal profongation algebra, it seems reasonable to suppose that 
finding new solutions is equivalent to evolving realizations of this algebra. We hope to 
encourage listeners to achieve a faithful realization of K2. 

III. Even More Complicated Vector-Field PDE's 

Quite an interesting system of vector-field-valued pde's has recently arisen in inves- 
tigations of a student, Denis Khetselius, who originally came to UNM from Dubna. His 
work on twisting, Petrov type N vacuum solutions lead him to the following (involutive) 
set of coupled equations: 

Cj£i - AiMj = [£i, Mj] , Vi, J = +, 0, - , 
£,+ = adb - edf , £_ = bda - fde , A+ = ud^ , A- = wdu , (3.1) 
jCo = [C+,C_], Ao = [A+,A-], 

where the vector fields Aij and, separately, the Si generate a realization of sl(2,C) 
in their (pseudopotential-type) variable spaces. Notice, of course, that the differential 
operators Cj and Ai also constitute realizations of the generators for sl(2,C). (The 
quantity s = af + eb is a characteristic for all the Ci. If one treats the 4 variables as 
complex, projective coordinates, in C"^, for the group manifold, s is the radius variable.) 

We can write the most general solution to this system of pde's, which is quite 
"messy." It can be done in a number of different, but equivalent, ways, depending 
upon the order of the integrations performed. However, one could hope for a much less 
coordinate-dependent approach to such a problem. To emphasize the meaning of this 
quest, we first study a slightly reduced version, obtained by assuming the /Aj to be 
constant: 

Cj£i = [Mj, £i\ = {adMj}£i . (3.2) 

These equations may be treated as saying that the £i are eigenvectors of the "total 
angular momentum" operators, with eigenvalue zero, taking the ad-action of the Aij as 
a realization of (the negative of) the usual 'spin'-operators for sl(2,C). Such a point of 
view ought, it seems, to generate "nice" expressions for the solutions; however, we have 
not yet seen them. 
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